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The NEMO-3 detector, which had been operating in the Modane Underground Laboratory from 
2003 to 2010, was designed to search for neutrinoless double [3 (Oiy/I/I) decay. We report final 
results of a search for Oz^yd^d decays with 6.914 kg of ^°°Mo using the entire NEMO-3 data set with 
a detector live time of 4.96 yr, which corresponds to an exposure of 34.3 kg-yr. We perform a 
detailed study of the expected background in the Oz/yd^d signal region and find no evidence of Oiy/d/I 
decays in the data. The level of observed background in the Oz/^d^d signal region [2.8 — 3.2] MeV is 
0.44 ± 0.13 counts/yr/kg, and no events are observed in the interval [3.2 — 10] MeV. We therefore 
derive a lower limit on the half-life of Oz^^d^d decays in ^°°Mo of Ti/ 2 ( 0 z/ydyd) > 1.1 x 10^^ yr at the 
90% Confidence Level, under the hypothesis of decay kinematics similar to that for light Majorana 
neutrino exchange. Depending on the model used for calculating nuclear matrix elements, the limit 
for the effective Majorana neutrino mass lies in the range {m^) < 0.33-0.62 eV. We also report 
constraints on other lepton-number violating mechanisms for Oz/yd/3 decays. 

PACS numbers: 23.40.-s, 21.10.-k, 27.60.+j 
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I. INTRODUCTION 


Since neutrinos are the only fermions that carry no 
electric charge, they can be represented by a Majorana 
field for which the distinction between matter and anti¬ 
matter vanishes. The Majorana nature of neutrinos could 
play a fundamental role in many extensions of the Stan¬ 
dard Model. For instance, the see-saw mechanism [T], 
which requires the existence of a Majorana neutrino, nat¬ 
urally explains the origin of small neutrino masses. A 
Majorana neutrino would provide a framework for lepton 
number violation, and in particular for the Leptogene- 
sis process [ 2 ] , which could explain the observed matter- 
antimatter asymmetry in the Universe. 

The observation of neutrinoless double [3 (Oz^/d/3) de¬ 
cay would prove that neutrinos are Majorana particles [3] 
and that lepton number is not conserved. The most com¬ 
monly studied mechanism of Oz//3/3 decay is the exchange 
of a Majorana neutrino. However, other mechanisms 
such as the existence of right-handed currents in the elec- 
troweak interaction, the exchange of supersymmetric par¬ 
ticles with R-parity violating couplings, or the additional 
emission of a Majoron particle, are possible. Except for 
the case of Majoron emission, the experimental signature 
of OuPP decays is the emission of two electrons with a to¬ 
tal energy Etot that is equal to the transition energy 
of the decay. 

For a given mechanism and isotope, the 0i'f3f3 decay 
half-life depends on the phase space factors and on the 
nuclear matrix element (NME). The decay half-lives of 
different isotopes can differ by a few orders of magnitude 
with large theoretical uncertainties of the NME calcula¬ 
tions. It is therefore essential to search for OuPP decays 
in several isotopes. 

The NEMO-3 detector [4] was operated from 2003 un¬ 
til 2010 in the Modane Underground Laboratory (LSM) 
to measure two-neutrino double p (2vpp) decays of seven 
isotopes in the form of thin foils and to search for Oupp 
decays. The full topology of double P decays is recon¬ 
structed by combining information from a calorimeter 
and a tracking detector that are both distinct from the 
source foils. We measure the contributions from different 
background processes separately by exploiting specific 
event topologies. The NEMO-3 design and its capacity 
to identify electrons, positrons, 7 rays, and a particles 
are unique in enabling us to reject background processes 
very efficiently. 

The isotope ^^^Mo represents the largest source sample 
in NEMO-3 with a mass of 6.914 kg and = 3034.40± 
0.17 keV [ 5 ]. A result based on a subset of the data 
had previously been published in [6]. We reported as 
rapid communications [7] the results of a search for OuPP 
decays for the entire data set, corresponding to a live time 
of 4.96 yr and an exposure of 34.3 kg-yr of ^^^Mo. In this 
Article, we describe this analysis in more detail. 

The NEMO-3 detector is introduced in Section [llj and 
the energy and timing calibration of the detector are de¬ 
scribed in Section III Selection criteria for Oupp candi¬ 


dates are given in Section |rVj The methodology and the 
results of the measurement of the different background 
components are presented in Section |V| Results of the 


search for OuPP decays are summarised in Section VI 


II. THE NEMO-3 DETECTOR 

The distinctive feature of the NEMO-3 detection 
method is a full reconstruction of the double P decay 
topology using tracking in three dimensions as well as 
calorimetric and timing information. It provides not only 
the total energy Etot of the two simultaneously emitted 
electrons, but also the single energy of each electron and 
their angular distribution at the emission point from the 
foil. A detailed description of the NEMO-3 detector can 
be found in [4]. 

The thin source foils with a density of 40-60 mg/cuP 
containing the active double P decay isotope are sur¬ 
rounded by a tracking detector comprising open drift cells 
and a calorimeter composed of plastic scintillators. The 
source foils are distributed over a cylindrical surface of 
about 20 m^, which is segmented into 20 sectors of equal 
size, as shown in Figure 



FIG. 1: A schematic view of the NEMO-3 detector, showing 
the double P source foils, the tracking chamber, the calorime¬ 
ter composed of scintillator blocks and PMTs, the magnetic 
coil and the shield. 


Several double P decay sources are installed in the de¬ 
tector. The main isotope used to search for OuPP de¬ 
cays is ^^^Mo with a total mass of 6.914 kg. Smaller 
amounts of other isotopes are mainly used to measure 
2iypp decays, comprising ^^Se (0.932 kg, 2 sectors), ^^^Cd 
(0.405 kg, 1 sector), ^^^Te (0.454 kg, 2 sectors), ^^^Nd 
(36.55 g), ^^Zr (9.4 g), and "^^Ca (7 g). In addition, 1.5 
sectors of natural tellurium, corresponding to 0.614 kg of 
Te 02 , and 1 sector equipped with pure copper (0.621 kg) 
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are used to perform measurements of backgrounds from 
processes other than double /3 decay. There are seven 
foil strips in each sector. The mean length of the strips 
is 2480 mm with a width of 65 mm for the five central 
strips and 63 mm for the two edge strips. 

There are two types of ^^^Mo foils, metallic and com¬ 
posite. The metallic foils were produced in vacuum by 
heating and rolling ^^^Mo mono-crystals in the form of 
foils. To produce the composite foils, thin and chemi¬ 
cally purified ^^^Mo powder was mixed with polyvinyl 
alcohol (PVA) glue and then deposited between mylar 
foils with a thickness of 19 fim. The metallic Mo foils 
were placed in Sectors 02, 03 and 04. There are also 
five additional strips in Sector 1 and two strips in Sec¬ 
tor 5. The total surface of metallic foils is 43924 cm^. 
The total mass of ^^^Mo in metallic foils is 2479 g and 
the average percentage of ^^^Mo enrichment is 97.7%. It 
corresponds to an average surface density of the metal¬ 
lic foils of 57.9 mg/cm^. The composite Mo foils were 
placed in Sectors 10, 11, 12, 13, 14, 15, and 16. There 
are also two additional strips in Sector 01 and three strips 
in Sector 05. The total surface of composite Mo foils is 
84410 cm^. The total mass of ^^^Mo in composite foils is 
4435 g and the average percentage of ^^^Mo enrichment 
is 96.5%. The total mass of components (Mo, PVA and 
Mylar) is 5569 g. It corresponds to an average surface 
density of the composite foils of 66.0 mg/cm^. 

On both sides of the source foils, a gaseous tracking 
detector comprising 6180 open drift cells operating in 
Geiger mode provides three-dimensional track informa¬ 
tion. We use a cylindrical coordinate system with the 
z-axis pointing upwards. The drift cells are oriented par¬ 
allel to the ^-axis and provide measurements of the trans¬ 
verse and longitudinal coordinates of the track. To min¬ 
imize multiple scattering, the gas is a mixture of 94.9% 
helium, 4% ethyl alcohol, 1% argon, and 0.1% water va¬ 
por for a total volume of about 28.5 m^. 

The basic cell consists of a central anode wire sur¬ 
rounded by eight ground wires. All the wires are 50 yum 
in diameter and 2.7 m long. On each end of the cell is 
the cathode ring. When a charged particle crosses a cell 
the ionized gas yields around six electrons per centime¬ 
ter. These secondary electrons drift towards the anode 
wire at a speed of around 1 to 2 cm/fis depending on the 
distance of the electrons to the anode. Measurements of 
these drift anodic times are used to reconstruct the trans¬ 
verse position of the particle in the cells. In the Geiger 
regime, the avalanche near the anode wire develops into 
a Geiger plasma which propagates along the wire at a 
speed of about 6.5 cm//is. The arrival of the plasma 
at the two ends of the wire is detected with the cathode 
rings mentioned above. The two propagation times of the 
plasma are used to determine the longitudinal position of 
the particle as it passes through the cell. 

To readout the drift cells, the analog Geiger signals 
from the anode wires and the two cathode rings signals 
are first amplified and then compared to anode and cath¬ 
ode thresholds. For signals exceeding the thresholds. 


the anode signal starts four Time-to-Digital Gonverter 
(TDC) scalers. The first three are for the anode and the 
two cathode contents which are measured with a 12-bit 
TDC and give times between 0 and 82 /is. The last TDC 
scaler is 17-bits, which can provide time measurements 
between 0 and 2.6 ms. It is used for delayed a particle 
tagging. The cathode TDCs are stopped by the cathode 
signals while the anode TDCs are stopped by a signal 
sent by the general trigger. 

The methode for the track reconstruction and its cali¬ 
bration, and the tracking performances are presented in 
[4]. The average transverse and longitudinal resolutions 
of the Geiger cells are 0.5 mm and 0.8 cm, respectively. 
If the two electron tracks from a double (3 decay are con¬ 
strained to originate from the same vertex in the foil, the 
transverse and longitudinal vertex resolutions, defined as 
the r.m.s. of the distance between the intersection points 
of the two individual tracks with the foil, are 0.6 cm and 
1.0 cm, respectively. These resolutions are sufficient to 
discriminate between decays from different source foils 
and isotopes. 

The energy and time-of-flight of particles are measured 
by polystyrene scintillators surrounding the tracking de¬ 
tector. We use the time-of-flight to discriminate between 
signal events emitted from the foil and background events 
where particles crossed the foil. To further increase ac¬ 
ceptance, the end caps (the top and bottom parts of the 
detector, named petals) are also equipped with scintil¬ 
lators in the spaces between the drift cell layers. The 
calorimeter is composed of 1940 optical modules, which 
consist of large scintillator blocks, with a typical size of 
20 X 20 X 10 cm^, coupled to low radioactive photomul¬ 
tipliers (PMTs). 

The analog PMTs signals are sent to both a low and 
a high threshold leading edge discriminator. If the PMT 
signal exceeds the lower level threshold it starts a TDG 
measurement and opens a charge integration gate for 80 
ns. The high threshold discriminator works as a one shot 
that delivers a calorimeter event signal to the trigger logic 
which reflects the number of channels that have exceeded 
the upper threshold. This level is used to trigger the 
system (first level trigger) if the desired multiplicity of 
active PMTs is achieved. The trigger logic then pro¬ 
duces a signal called STOP-PMT, which is sent to all 
the calorimeter electronic channels, to save their data. 
So the TDCs are stopped and the integrated charge is 
stored. Then digital conversions begin. At the same 
time, a signal is sent to the calorimeter acquisition pro¬ 
cessor, which permits the read out of the digitized times 
and charges for the active channels. The analog-to- 
digital conversions of the charge and the timing signal 
are made with two 12-bit ADCs. The energy resolution 
is 0.36 pC/channel (about 3 keV/channel) and the time 
resolution is 53 ps/channel. If any PMT signal exceeds 
the high level threshold then the TDC measurement and 
charge integration are aborted and the system resets after 
200 ns. 

The external wall of the calorimeter is equipped with 5- 
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inch PMTs and the internal wall with 3-inch PMTs, and 
the end caps with both PMT types. The avera ge energy 
resolution of the calorimeter cje/E = 5.8%/^^^(MeV) 
for the scintillat ors equip ped with 5-inch PMTs, and 
(JeIE = 7.2%/Y^E(MeV) for the scintillators equipped 
with 3-inch PMTs. 

Photons are identified as hits in the calorimeter where 
no electron track points at the scintillator block. The 
scintillator blocks with a thickness of 10 cm yield a high 
photon detection efficiency of 51% (33%) for a photon of 
energy 1 MeV (3 MeV) at a normal angle of incident. 

A solenoidal magnet surrounding the detector provides 
a magnetic field of 25 G used to discriminate between 
electrons and positrons with an efficiency of about 95% at 
an electron/positron energy of 1 MeV. An external shield 
with a thickness of 19 cm constructed of low-radioactive 
iron, a borated water shield, and a wood shield surround 
the detector to reduce background from external 7 rays 
and neutrons. Calibrations are performed by inserting 
calibrated radioactive sources into the detector through 
dedicated tubes installed between each sector in the plane 
of the foils. 

During the first data taking period, labeled Phase I, 
from February 2003 until October 2004, the dominant 
background to the 0uf3l3 signal was contamination from 
radon (^^^Rn) in the tracking chamber. Radon contami¬ 
nation in the tracking chamber is measured by detecting 
electrons from (3 decay of accompanied by a de¬ 

layed a particle from ^^^Po decay. To detect delayed 
a particles, every hit inside the wire chamber arriving 
with a delay of up to 700 //s is read out with dedicated 
electronics. The ^^^Rn activity of about 30 mBq/m^ in¬ 
side the tracking chamber during Phase I is caused by a 
low rate of diffusion of ^^^Rn from the laboratory hall, 
with an activity of around 15 Bq/m^, into the detector. 
This contamination was significantly reduced, by a fac¬ 
tor of about 6 , by the installation of a radon-tight tent 
enclosing the detector and a radon trapping facility in 
December 2004. The second data taking period between 
December 2004 until the end of running in December 
2010 (Phase II) therefore has a reduced radon gas con¬ 
tamination of around 5 mBq/m^. Data from both Phases 
are presented in this Article. 

The trigger conditions used for recording double (3 can¬ 
didate events require at least one PMT signal with an 
amplitude greater than 50 mV, corresponding to an en¬ 
ergy of > 150 keV deposited in the associated scintilla¬ 
tor, in coincidence with at least three hits in the tracking 
detector within a time window of 6 /is recorded in the 
same half-sector of the detector as the scintillator hit. 
Additional PMT signals with an amplitude of > 10 mV, 
corresponding to an energy deposit of > 30 keV, are also 
recorded if they coincide within a time window of 80 ns. 
The trigger rates of the data acquisition are about 7 Hz 
for Phase I and about 5 Hz for Phase-H. The dead time of 
the data acquisition is measured to be 1 % and is treated 
as an inefficiency. 

Monte Carlo (MC) simulations are performed with 


a GEANT3-based [8] detector simulation using the 
decayO [ 9 ] event generator. The time-dependent sta¬ 
tus and conditions of the detector and its performance 
are taken into account in the detector simulation. 

In this Article, we present a search for 0i'l3f3 decays 
using data recorded between February 2003 and October 
2010, with a live time of 1.02 yr in Phase I and 3.94 yr 
in Phase H, and a total mass of 6914 g of ^^^Mo in the 
form of metallic and composite foils. This corresponds 
to a total exposure of 34.3 kg-yr. 


III. CALIBRATION 
A. Energy scale calibration and resolution 

Absolute energy calibrations of the calorimeter optical 
modules were carried out every month using ^^^Bi sources 
which provide internal conversion electrons with energies 
of 482 keV and 976 keV from the K lines, with branching 
ratios of 1.5% and 7.1%, respectively. Each calibration 
run has a length of about 24 hours. In addition, a ded¬ 
icated long run was performed using a ^^Sr source since 
the end point of the (3 spectrum of a daughter nu¬ 
cleus of ^^Sr, provides an additional high-energy point at 
an energy of 2279 keV. 



E3 (MeV) 


FIG. 2: Energy spectrum of a typical scintillator block, mea¬ 
sured with the ^°^Bi calibration sources and summed over 
all the calibration runs. The data points are compared to a 
histogram of the energy spectrum calculated by the MC simu¬ 
lation. The peaks correspond to the energies of electrons from 
the main 482 keV, 976 keV and 1682 keV internal conversion 
K lines of ^°^Bi. 


The response of each scintillator block to electrons with 
an energy of 976 keV is measured as a function of the im¬ 
pact position of the electron track on the entrance surface 
of the scintillator using ^^^Bi calibration runs. A depen- 
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dence on impact position was previously observed with 
data obtained with the electron spectrometer during the 
NEMO-3 calorimeter assembly. The impact position is 
sampled by dividing the entrance surface of the scintilla¬ 
tor blocks in 3 X 3 equal squares for the blocks equipped 
with 3-inch PMTs, and 5x5 equal squares for the blocks 
equipped with 5-inch PMTs, corresponding to 3 x 3 and 
5x5 corrections points respectively. The impact correc¬ 
tions are small for the scintillator blocks equipped with 
3-inch PMTs, typically 1% —2%, but they can increase up 
to 10% for 5-inch PMT scintillator blocks. This effect is 
corrected offline by applying different impact correction 
factors for each scintillator block type. 

The linearity of the PMTs has been verified with a ded¬ 
icated light injection test during the construction phase. 
Upper limits on the non-linearity of the PMTs were found 
to be < 1% for energies < 4 MeV, corresponding to 
the energy range of interest for double (3 decay mea¬ 
surements. It is shown by Monte-Carlo simulations that 
a non-linearity lower than 1% has no effect in the final 
f3f30u analysis. 

The linear fit combining the energy calibration ob¬ 
tained with the two ^^^Bi energy peaks and the end point 
of the (3 spectrum does not intersect with the origin, 
because the scintillator response for electrons at low ener¬ 
gies (below the energy threshold of 200 keV) is non-linear. 
The extrapolated energy offset at a charge of Qadc = 0 is 
on average 33 ± 3 keV. It is determined after subtracting 
the electronic pedestal of the Analog-to-Digital Convert¬ 
ers (ADCs) used to read out the PMTs and accounting 
for an impact point correction. This offset is taken into 
account in the energy calculation. It is shown by Monte- 
Carlo simulations that the uncertainty on the energy off¬ 
set measurement is negligible for the final analysis. 

An example of a linear fit for one counter can be found 
in Ref [4]. 

The rare internal conversion electron K line of ^^^Bi 
with an energy of 1682 keV has a small branching ratio 
of 0.02%. It is used to determine the systematic uncer¬ 
tainty on the energy scale from the difference between 
the reconstructed peak position in data and MC simula¬ 
tion, which is < 0.2% for 99% of the optical modules. It 
is shown by Monte-Carlo simulations that a 0.2% uncer¬ 
tainty on the energy scale is negligible for the final PpOu 
analysis. The remaining optical modules of the calorime¬ 
ter with incorrect reconstruction of the energy peak are 
rejected in the analysis. A typical energy spectrum mea¬ 
sured with a single optical module is shown in Figure 

Figure shows the average energy resolution as a func¬ 
tion of running time for the different types of scintillator 
blocks and PMTs. The resolution at an electron energy 
of I MeV ranges from cjeIE = 5.7% to 8.0%, depending 
on the type of block and the data taking period. A de¬ 
terioration of the energy resolution of 0.03%-0.05% and 
0.06%-0.I4% per year is observed for the blocks equipped 
with 5 and 3-inch PMTs, respectively. This drift might 
be caused by the residual helium concentration in the 
air surrounding the detector that leads to after-pulsing 



FIG. 3: Average energy resolution geIE measured at F = 
I MeV for the different types of scintillator blocks and sizes 
of PMTs as a function of NEMO-3 running time. Here, IN 
refers to the calorimeter blocks with 3-inch PMTs located in 
the central tower (inner wall of the calorimeter), LI, L2 and 
L3 refer to 3-inch PMTs located on the upper and lower end 
caps, EC and EE refer to 5-inch PMTs located on the external 
wall, and L4 to 5-inch PMTs on the upper and lower end caps 
of the calorimeter (see Figure 3 in [1] for the exact location 
of the different types of scintillator). 


of the PMTs. The helium concentration in the central 
tower of the NEMO-3 detector, where most of the 3-inch 
PMTs are located, is higher than in other regions of the 
detector, which could explain the larger drift in this re¬ 
gion. 

The expected full width at half maximum (FWHM) of 
the spectrum of two electrons energy sum in decays 

is 350 keV. It is a convolution of the energy resolution 
of the calorimeter and of the non-Gaussian fluctuations 
in the electron energy loss, which occur mainly in the 
source foil and to a lesser extent in the tracking detector. 
In the absence of energy loss fluctuations in the foil, the 
expected FWHM would be about 250 keV. 

After close to eight years of stable operation of the 
experiment, fewer than 10% of PMTs had to be turned 
off because they displayed unstable gain or noisy signals. 
The fraction of dead PMTs as a function of the NEMO- 
3 running time is presented in Figure The fraction 
of PMTs with noisy signals in the recorded data is es¬ 
timated by measuring the random coincidence rate of 
scintillator hits with a constant timing distribution. The 
same fraction of PMTs is randomly rejected in the MC 
simulation, leading to a reduction of the 0ul3f3 detection 
efficiency of 0.9% in Phase I and 2.4% in Phase II. 
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FIG. 4: Fraction of dead PMTs (in %) as a function of NEMO- 
3 running time for 3 and 5-inch PMTs, and for all PMTs. 


B. Laser survey 

The stability of the PMT gains between two consecu¬ 
tive absolute ^^^Bi calibration runs is maintained using 
dedicated laser runs which were performed twice daily. 
The laser beam is split and transmitted to two different 
devices to calibrate the 3 and 5-inch PMTs separately. A 
description of the laser system is given in [4]. Data tak¬ 
ing is divided into successive laser survey periods that are 
separated by major incidents such as a general shutdown 
of the high voltage crates or any other event that could 
cause a discontinuity in the operating conditions of the 
PMTs. 

The laser survey measures the time dependence of the 
average response of all PMTs of the same size to monitor 
the variation of gains. The mean energies {es{t)) and 
(e 5 (t)) of the two different sets of 3 and 5-inch PMTs 
calculated for each laser run at a recording time t are 
given by 


('».=(•)) =- rkjif, -• m 

where the sum extends over the 3 or 5-inch PMTs. Here, 
QABc{k,t) and are the recorded charge after 

pedestal subtraction and the laser calibration constant 
for the PMT labeled k and for the laser run recorded 
at time t. The numbers of 3 and 5-inch PMTs recorded 
during a laser run are (t) and (t), respectively. The 
parameters r]{k,t) are calculated for each PMT 




9k'^''"{t)QADG{k,t) 

(e3,5(i)) 


( 2 ) 


depending on its type. The parameters r]{k, t) are divided 
by {r]o{k)), which is the mean value of rj{k,t) during the 


associated absolute energy ^^^Bi calibration run, to cal¬ 
culate the final laser correction factor of 

C{k,t) =ri{k,t)/{rio{k)). (3) 


(/> 



Laser correction factor 

FIG. 5: Distribution of the laser correction factors calculated 
for all laser runs and for all stable PMTs used in the analysis. 


The time dependence of the correction factors is anal¬ 
ysed to characterize the level of stability of the PMT 
gains during data taking. A large change of the correc¬ 
tion factor or discontinuities during a period between two 
absolute energy calibrations are interpreted as an insta¬ 
bility, and the corresponding PMT and associated events 
are rejected for that period. For each PMT, we esti¬ 
mate its stability during that period, by determining the 
number of laser runs for which the correction factor de¬ 
viates 5%. During the entire data taking period, 82% of 
PMTs are considered to be stable. Taking into account 
that more than 90% of data are recorded with a reliable 
laser survey, the efficiency to select a double (3 event is 
reduced by 25.2% in Phase II when the laser survey is 
applied, as reported in Table [TT| The efficiency reduction 
is 38.6% in Phase I because of a less stable laser during 
this first phase of data. The distribution of the laser cor¬ 
rection factors for all laser runs is shown in Figure for 
the stable PMTs. 

The reliability of the laser survey procedure is vali¬ 
dated by analysing a pure sample of electrons with an 
energy close to the end point of (3/3 = 3.27 MeV in the /3 
energy spectrum of ^^^Bi decays occurring in the track¬ 
ing chamber. Any excess in data over the MC expecta¬ 
tion around would be a sign of unstable PMT gains. 
The events are selected by requiring electrons in coinci¬ 
dence with a delayed a track from the ^^^Bi-^^^Po cas¬ 
cade (“BiPo events”)- 

The entire data set is used in this analysis. The se¬ 
lection of BiPo events is similar to the one used for 
the radon background measurement, described in Sec¬ 
tion VB| Here, only BiPo events with a vertex inside 
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the tracking chamber are selected, and the electron track 
length is restricted to > 45 cm. Electrons crossing the 
source foils are rejected, since they could have lost energy 
in the foils. 

To minimize the proportion of re-firing Geiger cells, we 
require the a delay time to be > 140 /is for events with 
only one delayed hit, and > 70 /rs for > 1 hits. The 
delay time distributions are analysed separately for each 
Geiger plane, on each side of the source foils, and for 1 , 
2, 3, and > 3 delayed hits. Regions within the tracking 
chamber with a significantly increased fraction of random 
coincidences or re-firing cells are excluded. These effects 
are therefore negligible in the selected data set. 



Phase 1 

Phase 2 


Data MC 

Data MC 

No laser correction 

8 2.32 ±0.32 

9 0.77 ±0.15 

With laser correction 

2 1.72 ±0.28 

1 0.50 ±0.11 


TABLE I: Numbers of BiPo events with Ee > 3.4 MeV 


The electron energy spectra obtained from this analy¬ 
sis are shown in Eigurej^ for Phases I and II separately, 
before and after applying the laser corrections. They are 
compared to the expected background from the MG sim¬ 
ulation, assuming the ^^^Bi activity on the surfaces of 
wires and foils, and the ^^^Bi activity inside the foils de¬ 
scribed in Section [Vl The number of data and MG events 
with Ee > 3.4 MeV, without laser correction and after 
applying the laser correction is given in Table [l| Spurious 
events observed beyond the end-point are well rejected 
after applying the laser survey. It demonstrates the re¬ 
liability of the laser survey to reject false high energy 
events with a wrong recorded energy. 


C. Timing calibration and time-of-flight 

Time-of-flight measurements are used to discriminate 
between two-electron events from double P decays emit¬ 
ted from the source foil and events where an external elec¬ 
tron crosses the detector and foil. The crossing electron 
in these events could be reconstructed as two separate 
tracks with a common vertex. 

The time calibration of the optical modules takes into 
account both the individual absolute time shift of each 
optical module and a time-vs-charge dependence induced 
by the effect of leading edge discriminators. The cali¬ 
brated time, t(i), used for a time-of-flight calculation for 
counter number i is: 

t{i) = tdc{i) - ts{i) - (4) 

where tdc{i) is the TDC measurement, ts(i, t) is the time 
shift, and is the time-charge correction function, 

which correct the measurement of the TDG as a function 
of the charge Q (with formula given in [4]). 


The absolute time shifts are measured individually for 
each optical module, using (e“, 7 ) events selected from 
the absolute energy calibration runs carried out with 
^^^Bi sources using the relation 


^ _ Ejv,(^C/ + + fiQij)) - f{Q{i)) + ts{j)) 

~ Ni 

(5) 

where = ttof{j) — ttofii) is the difference between 

the calculated time-of-flights ttof of the electron and the 
7 , = tdc{j) — tdc{i), and Ni is the number of 

selected events with the optical module i hitted 

by the electron. 

The time-vs-charge correction functions f(Q) are mea¬ 
sured for each of the seven scintillator block types by us¬ 
ing crossing-electron events from a dedicated run with 
an external Am-Be fast neutron source. East neutrons 
are thermalized mostly in the scintillators. Then 7 are 
created by the capture of the thermalized neutrons in 
the copper walls. If a 7 produce an electron by Comp¬ 
ton effect in the scintillator, this electron can escape 
and crosses the tracking chamber, producing a crossing- 
electron event. The correction functions are measured 
using the relation 


f(Q(i)) = Atdc^’^ - (ts(i) - + f{Q{j)) (6) 


where AtdE^^ = ttof{i) + ttof{j) is the calculated time- 
of-flight for the electron to cross the tracking detector 
from the optical module j to the optical module i. The 
values of f{Q{i)) are groupped according to the seven 
scintillator block types, and then used to produce the 
time-vs-charge f{Q) distribution that is then fitted with 
a formula using four parameters pk . 


f{Q) = Pi 


P2 

P3 X \/Q + P4 


(7) 


Since the absolute time shifts and the time-vs-charge 
correction functions are both used in the two calibration 
relations and an iterative procedure is required to 
determine them. Eirst the absolute time shifts are calcu¬ 
lated according to Equation with initial values of the 
time-vs-charge correction functions obtained with laser 
runs and initial values of the time shifts set to zero. Then 
the time-vs-charge correction functions are calculated us¬ 
ing Equation!^ These new corrections functions are then 
used to calculate the absolute time shifts, and so one. 
Successive iterations are performed until a convergence 
is obtained. 

The daily laser surveys are used to identify and correct 
any variation of the TDG response. This laser timing cor¬ 
rection is calculated separately for each optical module 
and laser survey run. 

The average timing resolution of a scintillator hit is 
about 250 ps for a I MeV electron. 

The time-of-flight analysis is based on a comparison 
between the measured and expected time differences of 
the two scintillator hits. The expected time-of-flight is 
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FIG. 6: Energy spectrum of electrons from [3 decay of measured using BiPo events inside the tracking detector, without 
laser survey (a,b) and after laser survey (c,d) for Phase I (a,c) and Phase II (b,d). The data are compared to a MC simulation. 
The excess of electrons observed at Ee > 3.4 MeV in data are caused by PMTs with unstable gains. They are rejected by the 
laser correction (see Table |^. 


calculated assuming two different hypotheses: the ex¬ 
ternal hypothesis corresponding to a crossing electron 
and the internal hypothesis corresponding to two elec¬ 
trons being emitted simultaneously from the same ver¬ 
tex on the foil in a double (3 decay. The time-of-fiight 
calculation also accounts for the length of the tracks 
and the energy loss in the tracking detector. To cor¬ 
rectly take into account uncertainties on the timing mea¬ 
surement, we calculate separate probabilities for internal 
two-electron events (Pint) and external crossing-electron 
events (Pext)- The distributions of the difference AT be¬ 
tween the measured and theoretical time differences of 
the two scintillator hits, calculated assuming the internal 
hypothesis, is shown in Figure for the full sample of 
two-electrons events selected using all criteria described 
in Section IV, except the requirement on the time-of- 
fiight. The Pint distribution shown in Figure is con¬ 
stant above Pint = 1 %, as expected for double p decays, 
while the peak at Pint < 1 % corresponds to crossing- 
electron events. Internal double [3 events emitted from 
the source are centred around AT = 0 ns, while crossing- 
electron events from external background sources have 
I AT I > 3 ns. The r.m.s. of the AT distribution for 
Pint > 1% is 490 ps. 


IV. SELECTION OF DOUBLE [3 DECAY 
EVENTS AND EFFICIENCY 


Candidate double P decay events are selected by re¬ 
quiring exactly two electron tracks. Events with more 
than two tracks are rejected. 


• Each track must be associated with a scintillator 
hit, and the extrapolated track must hit the front 
face of the scintillator block and not the lateral 
side of petal blocks. The associated scintillator hits 
must be isolated, i.e., no hits are found in neighbor¬ 
ing scintillator blocks, and only a single track can 
be associated with the scintillator block. Events 
with a 7 candidate, defined by a scintillator hit that 
is not associated to a track, are rejected. 

• The two electron tracks must originate from a com¬ 
mon vertex in the ^^^Mo source foil. We therefore 
require that the transverse and longitudinal com¬ 
ponents of the distance between their intersection 
points with the foil are less than 4 cm and 8 cm, 
respectively. 

• To reject background from ^^^Bi decays near the 
foil, the number of unassociated hits in the tracking 
detector close to the vertex should not exceed one. 
When the two tracks are on the same side of the foil, 
there must be no unassociated hit on the opposite 
side of the foil close to the vertex. 

• The energy of each electron as measured in the 
calorimeter must be > 200 keV. 

• The curvature of both tracks must be negative to 
reject positrons. 

• The time-of-flight must correspond to the two elec¬ 
trons being emitted from the same vertex in the 
source foil, requiring Pint > 1 % and Pext < !%• 
To ensure a reliable time-of-flight measurement, 
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AT (ns) 

(a) 



(b) 


FIG. 7: Distributions of the difference AT between the mea¬ 
sured and expected time differences of scintillator hits for the 
internal hypothesis (a) and the internal probability Pint (b) 
for two-electron events. The superimposed shaded histogram 
shows events with Pint >1%. 


the track length of each track must exceed 50 cm. 
Events with hits in scintillator blocks from the in¬ 
nermost circle of petals are rejected. 


• Events with delayed tracker hits close to the elec¬ 
tron tracks are rejected t o red uce ^^^Bi and radon 
background (see Section VB). The delay time of 
these hits is required to be greater than 100 /is for 
events with only one delayed hit, and 40 /is, 20 /is, 
and 4 fis for events with 2, 3, or > 3 delayed hits, 
respectively. These criteria reduce the sensitivity 
to spurious hits in cells close to the electron track. 


• Events are rejected if a scintillator hit is linked to 
a PMT that has been flagged by the laser survey 
as having unstable gain. 

A typical double P event is shown in Eigure i Only 
events with an energy sum Ptot > 2 MeV for the two elec¬ 
trons are considered in the 0z/y5/3 search. The efficiencies 
to select OuPP events are calculated using the MC sim¬ 
ulation, and are given in Table [ll] after each successive 
application of the selection criteria. The 0ul3f3 signal se¬ 
lection efficiency is 11.3% for Phase I and II combined 


and Ptot > 2 MeV. It reduces to 4.7% in the energy 
window Ptot = [ 2.8 — 3.2] MeV. This reduction is due to 
the fact that the Ptot energy spectrum of the /3/30z/ signal 
peaks around 2.8 MeV, i.e. 200 keV below the theoretical 
( 5 / 3/3 value, because of the energy losses of the electrons in 
the foil and in the wire chamber. If the inefficiency due to 
noisy Geiger cells and unstable or dead PMTs is removed, 
these efficiencies increase to 20.3% for Ptot > 2 MeV, and 
8.5% in the energy window Ptot = [ 2.8 — 3.2] MeV. 

The uncertainty on the signal efficiency is determined 
using dedicated runs with two calibrated ^^^Bi sources, 
with a low activity of around 180 Bq, at four different lo¬ 
cations inside the detector. The runs were taken in March 
2004, June 2004 and April 2006. The two conversion elec¬ 
trons emitted simultaneously by the ^^^Bi sources are se¬ 
lected. The criteria to select the two electrons events are 
the same as the ones used to select the double (3 events, 
except that the energy of the electrons must correspond 
to the expected energy of the conversion electrons and the 
common vertex of the two electron tracks must originate 
from the calibration sources. The reconstructed ^^^Bi ac¬ 
tivities are in agreement with the nominal values within 
5%, which is consistent within the expected systematic 
uncertainty. 


Selection Griteria 

Ideal 

Phase I 

Phase II 

Trigger 

0.973 

0.973 

0.973 

Two tracks reconstructed 

0.480 

0.478 

0.462 

Track-scintillator association 

0.352 

0.348 

0.331 

Associated PMTs not dead 

0.352 

0.321 

0.288 

No extra scintillator hit 

0.313 

0.287 

0.258 

Scintillator correctly calibrated 

0.313 

0.281 

0.245 

Gommon track vertex in foil 

0.280 

0.251 

0.218 

Tracks have hits near foil 

0.273 

0.244 

0.211 

No extra prompt hits near vertex 

0.271 

0.242 

0.209 

Track length > 50 cm 

0.252 

0.225 

0.194 

Scintillator energy > 200 keV 

0.245 

0.219 

0.189 

Negative track curvature 

0.223 

0.199 

0.172 

Isolated scintillator blocks 

0.219 

0.195 

0.169 

No scintillator at petals near foil 

0.209 

0.186 

0.161 

Timing requirement 

0.206 

0.184 

0.159 

Reject a particles 

0.206 

0.184 

0.159 

Energy laser survey 

0.206 

0.113 

0.119 

Ptot > 2 MeV 

Ptot > 2.8 MeV 

0.204 

0.085 

0.111 

0.047 

0.117 

0.049 


TABLE II: Evolution of the efficiency as a function of 

the successive criteria of selection for Phase I and II. “Ideal” 
refers to the detector without any noisy Geiger cell neither 
unstable or noisy PMTs. 


V. BACKGROUND MEASUREMENTS 


The NEMO-3 detector is unique in its ability to iden¬ 
tify electrons, positrons, 7 rays and delayed a parti¬ 
cles by combining information from the tracking detec¬ 
tor, calorimeter, and the track curvature in the mag- 
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FIG. 8: Transverse and longitudinal view of a reconstructed double /3 data event. Tracks are reconstructed from a single 
vertex in the source foil, with an electron-like curvature in the magnetic field, and are each associated to an energy deposit in 
a calorimeter block. 


netic field. This allows the separation of different non 
double P background processes by exploiting differences 
in their event topologies and final states. We distin¬ 
guish three background components, as illustrated in 
Figure 1^ namely the external background, the internal 
background, and the background from radon. We first 
measure the external background. Then, the radon and 
thoron backgrounds inside the tracking detector are mea¬ 
sured, setting the external backgrounds to their measured 
values. Finally, the internal ^^^Tl and ^^^Bi contami¬ 
nations inside the PP source foils are determined, with 
all other backgrounds fixed. A full description of the 
background analysis and preliminary background mea¬ 
surements with part of the NEMO-3 data set are given 
in Ref. [10]. Here, we report the results of the background 
measurements using the full data set. 


A. External background 

External background is produced by the interaction of 
external 7 rays originating from the natural radioactivity 
of the detector outside of the source, by external neu¬ 
trons undergoing neutron capture that results in emis¬ 
sion of 7 rays, or by cosmic rays. If an external 7 ray 
is not detected by a scintillator, it can reach the source 
foil without being tagged. It can then mimic a /3/3 event 
by creating an e+e“ pair, if the two photons from a sub¬ 
sequent positron annihilation remain undetected or the 
sign of the positron track curvature is incorrectly mea¬ 
sured. Double or single Compton scattering followed by 
Mqller scattering can also lead to a double /5-like topol¬ 
ogy. The different mechanisms are illustrated in Eigurej^ 


We measure the external background using both exter¬ 
nal ( 7 , e~) and crossing-electron events, as illustrated in 
Eigure 10 External events are selected requiring 


one isolated scintillator hit, assumed to be from the 7 
ray, and one electron track coming from the source foil 
and associated with a different scintillator. The time dif¬ 
ference between the scintillator hits must agree with the 
hypothesis that an external 7 ray has hit the first scintil¬ 
lator block before producing a Compton electron in the 
foil. 


Crossing electrons leave a track that traverses the de¬ 
tector and is associated with a scintillator hit on either 
side with a time-of-flight and a curvature consistent with 
a crossing electron. In this topology, an external 7 hits 
the first scintillator block from outside and then creates 
an electron by Compton scattering in the last few mil¬ 
limeters of the scintillator closest to the tracking detector. 
This Compton electron crosses the detector including the 
foil before hitting the second scintillator, depositing its 
entire energy. 

The external background is modelled by fitting the 
data in both channels assuming contaminations of ^^^Bi 
from 238 u 208 rp| from decays, inside the 

PMTs, scintillators, iron shield and iron structure, cos- 
mogenic ^^Co inside the mechanical structure, and exter¬ 
nal 7 rays from the laboratory environment. 

The 20^T1 and ^^^Bi contaminations inside the PMTs 
are the dominant components of the external background 
in the range Etot > 2 MeV. Their activities have been set 
to the values quoted in our previous background measure¬ 
ment with part of the NEMO-3 data set [ 10 ]. Activities 
of other components in the MC simulation are fitted to 
the data using a combined fit to the distributions of the 










11 


Two electrons produced by an external y 



Two electrons produced by internal or ^°®T1 
contaminations inside the PP source foil 




P-decay + 
internal conversion 


P-decay + 
Moller 



P-decay + 
Compton 


Radon background 



FIG. 9: Schematic view of the different components of the two-electron background: the external background produced by an 
external 7 ray, the internal background produced by internal ^^^Bi and ^°®T 1 contaminations in the ^°°Mo source foil, and the 
radon contamination inside the tracking detector. 


electron energy E ^-, the 7 energy the sum of the en¬ 
ergy Eq- + E^^ and the angle between the reconstructed 
7 direction and electron track. 


Figure 11 shows the energy spectra of the Compton 
electrons for external ( 7 , e~) events and the energy mea¬ 
sured in the last scintillator block hit cross¬ 

ing electrons. The fitted MC background model agrees 
with the data and lies within the 10 % systematic un¬ 


certainty of the previous results obtained with a smaller 
data set HO]. It is also consistent with the radioactivity 
measurements of the detector materials using high-purity 
germanium (HPGe) detectors before installation [TO]. 

The neutron contribution to the external background 
is measured with dedicated runs performed with an Am- 
Be neutron source located outside of the shield. The 
data provide the energy spectra of Compton electrons 
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external (y ,e“) 
Y/ 



crossing e“ 
Y/ 



FIG. 10: The two events topologies used to measure the external background: external ( 7,0 ) events and crossing-electrons 
events. 


created by external neutrons in the ( 7 , e~) and crossing- 
electron channels. These spectra are then used in the 
fit of the external background model in Figure [m The 
contribution of neutrons to the external background is 
negligible for F^tot <2.6 MeV, which corresponds to the 
energy of the 7 line of ^^^Tl, but neutrons dominate at 
higher energies. The good agreement between data and 
expected background from neutrons shows that the mea¬ 
surement performed with the Am-Be source correctly 
emulates the expected external background induced by 
neutrons for E’tot > 2.6 MeV, and can be used to esti¬ 
mate the expected background in the energy range. 
Only six double P-\ike events with vertices in the ^^^Mo 
foils and 2.8 < F^tot < 3.2 MeV are observed in the 
Am-Be neutron data. With the normalization factor ob¬ 
tained from the fit of the external background in Fig¬ 
ure 11 , we obtain a negligible expected background rate 


of 0.03±0.01 events for the combined Phase I and II data 
sets in the energy range 2.8 < F^tot <3.2 MeV consistent 
with a 01/13f3 signal. The expected number of double /3- 
like events for E’tot > 4 MeV is 0.14 ± 0.03 and is also 
negligible. 

The neutron background model is further studied us¬ 
ing events with e+e“ pairs, where external neutrons 
are the only expected component of the background for 
F^tot > 4 MeV. The criteria to select e+e“ events are the 
same as the ones used to select two-electrons events (see 
Section IV), except that the curvatures of the two tracks 
are required to be of opposite sign. For E’tot > 4 MeV, 
we observe 2 e+e“ events, in agreement with the expec¬ 
tation of 1.1 ± 0.1 neutron events. The F^tot distribution 
for these events is shown in Figure [T^ 


B. Radon and thoron contaminations 

Radon and thoron are both found inside the track¬ 
ing detector. Radon (^^^Rn) with a half-life of T 1/2 = 


3.824 days and thoron (^^^Rn) with T 1/2 = 55.6 s are 
(a-decay isotopes that have ^^^Bi and ^^^Tl as daugh¬ 
ter isotopes in their respective decay chains. Radon and 
thoron emanate from the rock into the air, from where 
they diffuse into the detector and contaminate the in¬ 
terior of the tracking chamber. They can also emanate 
directly from the detector materials inside the tracking 
chamber. Subsequent a decays of these rare gases pro¬ 
duce ^^^Pb or ^^^Pb ions, which drift mainly to the cath¬ 
ode wires. If they are deposited on wires close to source 
foils, their decays can mimic di (3(3 decay, as illustrated in 
FigureContamination from thoron is much lower than 
from radon since the shorter half-life makes it less likely 
for thoron to emanate and diffuse into the detector. 

The radon contamination is measured by detecting 
BiPo events, where the electron from p decay of ^^^Bi, 
a daughter of ^^^Rn, is followed by a delayed a parti¬ 
cle from the decay of ^^^Po, which has a short half-life 
of 164 /is. Additional photons may also be emitted and 
detected. A BiPo event in the NEMO-3 detector is identi¬ 
fied by requiring an electron track inside the wire cham¬ 
ber associated with a scintillator hit, and at least one 
delayed hit in the tracking chamber close to the emission 
point of the electron, due to the delayed a particle. The 
delay time is required to be at least 100 /is for events 
with only one delayed hit, and at least 40, 20, and 4 /is 
for events with 2, 3 and > 3 delayed hits, respectively, 
to reject hits where electrons have caused neighboring 
Geiger cells to re-fire. Applying these criteria, the mean 
efficiency to select a BiPo event produced on the surface 
of a wire is estimated by MC simulations to be 23%. 

The time distribution of delayed tracks, shown in Fig¬ 
ure is used to demonstrate the purity of the event 
selection. We fit the sum of an exponential function 
and a constant term accounting for random coincidences 
to the data distributions, assuming a ^^^Po half-life of 
Ti /2 = 164 /is. For Phase II the fits are applied to delay 
times larger than 140 /is for events with only one delayed 
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FIG. 11: Result of the fit of the external background to data for the total ^°°Mo exposure of 34.3 kg-yr, for the electron energy 
Ee of external ( 7 , 0 “) events (a,c) and the energy measured in the last scintillator block hit in crossing-electron events 
(b,d). The distributions are shown separately for Phases I (a,b) and II (c,d). SC K40 corresponds to impurities inside the 
scintillators. Lower panels show residuals between data and expected background, normalized to the Poisson error, ignoring 
bins with 0 events. 


hit in the tracking detector, and 80 /is and 60 fis for 
events with 2 or > 2 delayed hits, respectively. Slightly 
lower minimum delay times are used for Phase I. The very 
small excess of events over the extrapolated curve at low 
delay time provides the fraction of re-firing Geiger cells, 
and the constant term provides the fraction of random 
coincidences. The contribution of random coincidences 
and Geiger re-firings, given in Table m depend on the 
number of delayed hits and the data taking period. In 
all cases, they are found to be negligible. 

This method allows a daily measurement of the radon 
activity inside the tracking detector. The average radon 
activity is about 30 mBq/m^ in Phase I and about 
5 mBq/m^ in Phase IT Figure [l4| shows the spatial dis¬ 
tribution of vertices for BiPo events that either originate 
on the foils or on one of the first two layers of Geiger 
cells inside the tracking chamber. The activity is larger 


in Sector 03, which hosts a ^^^Mo source, than in other 
sectors. The radon model used for the background sim¬ 
ulation includes the contributions of ^^^Bi deposited on 
the surface of wires and on the surface of foils. 

The systematic uncertainty on the ^^^Bi background 
contribution caused by radon contamination is domi¬ 
nated by the uncertainty on the efficiency of the track¬ 
ing chamber to detect a delayed a decay of ^^^Bi. It 
is estimated by independently measuring the activities 
of the isotope ^^^Bi using (e“,a) and (e“, 7 ) events. A 
large fraction of the ^^^Bi P decays are accompanied by 
a high energy 7 ray emitted from the same point inside 
the tracking chamber. These events are contam¬ 

inated both by external 7 rays that Compton scatter on 
the wires of the Geiger cells, and by (/3, 7 ) emitters in the 
wires. To suppress this background, only events with a 
7 energy > 1 MeV are selected. 
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'°Mo, eV, Phase I + Phase II 
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FIG. 12: Distribution of F^tot for e^e“ pair events consis¬ 
tent with being emitted from ^°°Mo foils for the entire data 
set. The data are compared to the sum of the expected back¬ 
ground from external neutrons, 2z//3/3 events, and the other 
background components. 


Number of Delayed Hits 

1 

> 1 

Random Goincidences 

Phase I 

< 0.03% < 2.7% 

Refiring 

< 0.5% 

< 2.6)% 

Random Goincidences 

Phase II 

< 0.05% (I.I ±0.3)% 

Refiring 

< 0.7% 

< 0.7% 


TABLE III: Contribution of random coincidences and Geiger 
refirings in the selection of BiPo events used for the Radon 
measurement, for the high radon period (Phase I) and the low 
radon period (Phase II), requiring either exactly one or several 
delayed Geiger hits. Upper limits are given at 90% C.L. 


in the case of ^^^Tl with a conversion rate of 0 . 2 % for 
the 2615 keV 7 ray, which produces a conversion elec¬ 
tron with an energy of 2527 keV. Other processes are 
M 0 ller scattering of the /3-decay electrons in the source 
foil, or p decay to an excited state followed by a 7 under¬ 
going Compton scattering, which can be reconstructed 
as two-electron events if the 7 is not detected. 


The ^^^Bi measurement using (e“, 7 ) events suf¬ 
fers from larger background and has an approximately 
three times smaller detection efficiency compared to the 
method using delayed tracks. It is sensitive to the sys¬ 
tematic uncertainties on 7 detection, but it is not affected 
by systematic uncertainties on the a detection efficiency. 
The ^^^Bi and radon measurement using (e“, 7 ) events 
agree within 10 % with the result using an electron and a 
delayed a track [TO]. 

The 2^^T1 activity from thoron inside the track¬ 
ing chamber is measured using (e“, 77 ) and (e“, 777 ) 
events (see next Section). The ^^^Tl activity is about 
0.1 mBq/m^, both in Phase I and in Phase II. Taking 
into account the branching ratio of 36% for producing 
208Ti lYi the 232decay chain yields a thoron activ¬ 
ity of about 0.3 mBq/m^. The MC simulations predict 
that this thoron activity leads to a background for two- 
electron events with F/tot > 2 MeV that is a factor of 50 
smaller than the background originating from radon for 
Phase I, and a factor of 8 for Phase II. The 208xi contri¬ 
bution is therefore negligible in the 0ul3f3 energy region, 
and for decays with F/tot >2.8 MeV. 


C. Internal backgrounds 

Internal backgrounds originating from radioactive con¬ 
taminants inside the source foils are mainly due to [3 
decay of ^^^Bi with = 3.27 MeV and ^osrp^ with 
Qf 3 = 4.99 MeV. The two isotopes are products of the 
and 232xji decay chains, respectively. As illustrated 
in Figure the presence of ^^^Bi and 208xi can mimic 
/3/3 events by a /3 decay accompanied by an internal con¬ 
version electron process. This is the dominant channel 


1 208 contamination in the source foils 


The (3 decay of 208 x 1 is usually accompanied by two or 
three 7 rays. The 208 x 1 contamination inside the sources 
foils is therefore measured by selecting internal (e“, 77 ) 
and (e“, 777 ) events defined as one electron track origi¬ 
nating from the source foil that is associated with a scin¬ 
tillator hit, and two or three isolated scintillator hits. The 
time-of-flight must be consistent with the hypothesis that 
all particles are emitted from the track intersection with 
the foil. 

We require that the energy of the electron is in the 
range 0.2 < E^- < 1.5 MeV, > 0.2 MeV for all 
7 energies, and that the sum ^E^ < 3.5 MeV. The 
condition 

E^- (MeV) > (4(MeV) - 1.5 x X] ^7 (MeV)) ( 8 ) 


rejects ^^^Bi background. The highest energy photon 
must have E^ > 1.7 MeV to select the 2615 keV 7 line. 
Finally, we require Pint > 0.05, Pext < 0.01, and the z 
coordinate of the emission vertex of the electron must sat¬ 
isfy l^l < 120 cm. The distributions of Pg-? 
the total energy Pg- E^E^ are shown in Figure 15 The 
thoron and radon activities inside the tracking chamber 
are set to the values obtained from the prior measure¬ 
ments described in Section IV Bl 

The measured 208 x 1 activities of the metallic and com¬ 
posite Mo source foils, and of the copper and tellurium 
foils are given in Table [rV| for both event topologies com¬ 
bined. The data are in agreement with the upper limits 
from the HPGe measurements of 208 x 1 activities, prior 
to the installation of the foils in the detector. The two 
event topologies, ( 6 ^, 77 ) and (e“, 777 ), give consistent 
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NEMO-3 Phase I 

single delayed hit multipie deiayed hits 




NEMO-3 Phase II 

single delayed hit multipie delayed hits 




FIG. 13: Time distribution of delayed a tracks, measured for BiPo decays emitted inside the tracking detector, for Phase I ((a) 
and (b)) and Phase II ((c) and (d)), and for single delayed Geiger hit ((a) and (c)) or multiple delayed Geiger hits ((b) and 
(d)). The distributions are fitted by the sum of an exponential function with Ti /2 set to the ^^^Po half-life of T 1/2 = 164 /xs 
and a constant term accounting for random coincidences. 


results when analysed separately. The ^^^Tl activities of 
the copper and tellurium foils are used in section fV^ for 
the validation of the background model. 


The systematic uncertainty on the ^^^Tl activity is de¬ 
termined by using two radioactive sources (the iso¬ 
tope U is a parent of ^^^Tl). The 208xi activities of the 
sources are first calibrated by gamma spectroscopy with 
a coaxial HPGe detector, by measuring the intensity of 
the 7 line emitted in the decay of ^^^Pb to ^^^Bi with an 
energy of 238 keV, while the two 7 lines emitted in the 
decay of 20871 ^ with energies of 583 keV and 2615 keV 
are used to check the results. The HPGe detection effi¬ 
ciency is determined with a calibrated source that 

has an activity known to within 0.5%, and using a MC 
simulation of the setup. The sources are measured at 
four different distances between the source and the Ge 
crystal. The four activities obtained for each distance 
are combined to obtain a total statistical uncertainties 
of 0.7% and a systematic uncertainty of 3%. The two 


calibrated 232 u gources are then temporarily introduced 
into the NEMO-3 detector through the calibration tubes. 
We select (e“, 77 ) and ( 6 ^, 777 ) events and fit the ac¬ 
tivities of the two sources using a MC simulation of 232u 
decays. The results are given in Table |V| The largest 
sources of systematic uncertainty are the knowledge of 
the exact location of the sources (3%) and the kinematic 
selection criteria ( 6 %). This systematic uncertainty is 
estimated by allowing a variation of the energy require¬ 
ments, considering tracks that traverse only a single sec¬ 
tor, tracks only on the inner or outer side of the foils, 
and by accepting or rejecting scintillator blocks with an 
energy <150 keV. 

The results of the in-situ NEMO-3 and the HPGe mea¬ 
surements shown in Table lV| are consistent within their 
systematic uncertainties. We assign a systematic uncer¬ 
tainty of 10% to the 208 rpi activity measurement, corre¬ 
sponding to the larger difference between the in-situ and 
the HPGe measurements obtained for the second 232u 
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FIG. 14: The spatial distribution (vertical coordinate versus sector number) of the emission vertex of detected 
decay cascade events emitted inside the tracking detector close to the source foils for Phase IL Left (right) correspond to events 
with an emission vertex on the internal (external) side of the source foil, (a) and (b) correspond to a vertex on the foil or on 
the wires of the first layer of Geiger cells close to the foil, (c) and (d) correspond to a vertex on the wires of the second layer of 
Geiger cells, and (e) and (f) correspond to a vertex on the wires of the third layer of Geiger cells. The external side of Sector 
13 is not represented because of noise observed for Geiger cells in this zone. 


Source Foil 


iVobs Nb S/B 


e A A (HPGe) 
(%) (/xBq/kg) (/xBq/kg) 
(90% C.L.) 


'^^Mo Metal. 

823 

281 

1.93 

2.05 

87 

± 

4 

< 

100 

^°°Mo Gomp. 

2241 

617 

2.63 

2.15 

12 $ 


3 

< 

170 

Gopper 

75 

60 

0.25 

1.82 

11 

± 

3 

< 

33 

130rp^ 

563 

155 

2.64 

2.54 

206 

± 

10 

< 

500 

Te-nat 

741 

121 

5.14 

2.18 

301 

± 

12 

< 

830 


TABLE IV: Numbers of observed (e“, 77 ) and ( 6 ^, 777 ) 
events (Nobs), expected number of background events (Nb), 
signal-to-background ratio, ^°®T 1 signal efficiency (e), and 
measured ^^^Tl activity of the ^°°Mo metallic (Metal.) and 
composite (Gomp.) foils, the copper, ^^°Te and natural Te 
foils. The activities of the foils are compared to the HPGe 
measurements performed before their installation. Only sta¬ 
tistical uncertainties are given. 


source. 

The 2^^T1 background measurement is validated by us¬ 
ing the two-electron channel with at least one associated 
7 ray emitted in time from the source foil 



Activity (Bq) 

Source ( 1 ) Source ( 2 ) 

NEMO-3 

HPGe 

7.36 ± 0.03 ± 0.52 14.56 ± 0.05 ± 1.02 
7.79 ± 0.04 ± 0.21 15.91 ± 0.09 ± 0.43 


TABLE V: The ^^^Tl activities from sources obtained 
with the NEMO-3 detector and with HPGe 7 spectrometers. 


In the region where the sum of the two electrons energies 
satisfies F^tot >2.6 MeV, ^^^Tl contamination inside the 
foil dominates, whereas 2uP/3 decays are strongly sup¬ 
pressed by the selection criteria. Figure shows the 
total energy of two electrons F^tot for events 

for the entire ^^^Mo data set. The normalisations of the 
different background components are set to the previ¬ 
ously measured values and are not fitted to this distribu¬ 
tion. The data are in good agreement with the expected 
background, which is dominated by ^^^Tl contamination 
inside the foils. We observe 7 events in the ^^^Mo foils 
in the interval [2.8 — 3.2] MeV whereas 8.8 events are 
expected from the simulation. This independent check 
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NEMO-3 composite foils, e2Y & eSy 



E, (MeV) L[E^] (MeV) E^+LLE^] (MeV) 


NEMO-3 ^°°Mo metailic foiis, e2y & e3y 



E 3 (MeV) Z[E^] (MeV) (MeV) 

FIG. 15: Distributions of the energy of the electron, E^- , the energy sum and E^- + X) ^7 using (e“, 77 ) and (e“, 777 ) 

events for the combined ^°°Mo data set. The top panels show the composite and the bottom panels the metallic foils. The data 
are compared to the sum of the expected background from MC simulations and the htted ^^^Tl activity inside the ^°°Mo foils. 


validates the estimation of the ^^^Tl activity inside the 
foils within relatively large statistical uncertainties. 


2. contamination in the source foil 

The ^^^Bi contamination inside the source foils is mea¬ 
sured by analysing the distribution of the length of the 
delayed a tracks in BiPo events. It allows the discrim¬ 
ination of the ^^^Bi contamination inside the foils, and 
inside the mylar for composite foils, from the dominant 
radon background close to or on the surface of a foil. 

The criteria for the selection of the BiPo events are 
similar to the selection used for the radon activity mea¬ 
surement, except that the common vertex of the electron 
track and the delayed a track must be in the foil or in 
the first layer of wires of the tracking chamber. The ^^^Bi 
contamination inside the source foils is found by fitting 
the distribution of the delayed a track length, taking into 
account the other unknown activities as free parameters 


in the fit. These parameters are the ^^^Bi activities from 
radon deposition on the surface of the source foils and 
on the surface of the two closest layers of wires. Only 
Phase II data are used to reduce the radon background. 


The results of the fit are shown in Figures. EZl and 
for the ^^^Mo composite and metallic foils, respectively. 
The resu lts of the ^^^Bi activity measurement are given 
in Table VI for ^^^Mo foils, and also for copper, ^^^Te, 
and natural tellurium foils. They are in agreement with 
the upper limits obtained from HPGe measurements. 

The measured ^^^Bi contamination is checked by se¬ 
lecting two-electron events emitted from the ^^^Mo foils, 
where an associated delayed a track is emitted from the 
two-electron vertex (e“e“, a). This channel is dominated 
by radon background close to the foil and by ^^^Bi con¬ 
tamination from inside the foil. The criteria to select 
the two electrons are the same as those used for the se¬ 
lection of double /3 decay events (see Section IV). The 
criteria to select the delayed a track are identical to 
those used for the radon background measurement. Us- 
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FIG. 16: Distribution of the total energy of two electrons F^tot 
in the {e~e~,N^) channel for the ^°°Mo data set compared 
to the expected background from ^°®T1 contamination inside 
the foils and to the total expected background. The nor¬ 
malisations of the different background components are not 
htted, but set to the measured values. No event is observed 
for F^tot >3.7 MeV. 


Source Foil 

Activity 

Foil 

(mBq/kg) 

Activity 

Mylar 

(mBq/kg) 

A (HPGe) 
Foil+Mylar 
(mBq/kg) 

A (HPGe) 
Mylar 
(mBq/kg) 

^^^Mo Gomp. 

0.31 ±0.04 

1.05 ±0.06 

< 0.34 

< 0.67 

^°°Mo Metal. 

0.06 ±0.02 

No mylar 

< 0.39 

No mylar 

Gopper 

0.16 ±0.04 

No mylar 

< 0.12 

No mylar 

ISOrpg 

0.41 ±0.06 

1.81 ±0.17 

< 0.67 

3.3 ±0.5 

Te-nat 

0.37 ±0.05 

1.11 ±0.17 

< 0.17 

1.7 ±0.5 


TABLE VI: Measured ^^^Bi activity of the ^°°Mo metallic, 
^°°Mo composite, copper, ^^°Te, and natural Te source foils, 
compared to the HPGe measurements performed before their 
installation. Only statistical uncertainties are given. The 
fraction of the mylar mass relative to the total mass of the 
foil is in the range 5%-10%, depending on the foil. 


ing all foils, we observe six events with a {e~e~, a) 

topology in the energy range for the two electrons of 
^tot = [ 2.8 — 3.2] MeV in the combined Phase I and 
II data, while 9.4 ± 0.4 events are expected from simu¬ 
lations. Within large statistical uncertainties, this result 
confirms the prediction for the ^^^Bi background contri¬ 
bution in the 0i'f3l3 signal region. 


D. Validation of background model with copper 
and tellurium foils 


The complete background model is validated by select¬ 
ing two-electron events emitted from the copper, natural 
tellurium, and ^^^Te foils (Q/ 3/3 = 2527.518 ± 013 keV) 
using the criteria described in Section IV The data cor¬ 


respond to an exposure of 13.5 kg-yr. The internal con¬ 
taminations of these foils in ^^^Tl and ^^^Bi are mea¬ 
sured using the same methods as those used for the Mo 
foils (see sectiomVClmd VC 2 ). Results of the inter¬ 


nal contaminations measurements are given in Tables IV 


and[VI| Figure [1^ shows the distributions of the sum of 
the energies of the two electrons for E’tot > 2 MeV, and 
Table VII gives the number of events with F^tot > 2 MeV. 
The observed numbers of two-electron events agree with 
the expectation from the MC simulation calculated us¬ 
ing the background model, which is dominated by radon 
background. The number of 2i'/3/3 decays of ^^^Te in this 
energy region is expected to be negligible m- In the 
full data set, only 3 events with two electrons from the 
sectors containing copper, ^^^Te, and natural tellurium 
foils remain in the energy region F^tot = [ 2.8 — 3.2] MeV, 
compared to a MC expectation of 3.6 ± 0.2 events. 


Data Set 

Phase I 

Phase II 

Gombined 

External Background 

4.77 ±0.48 

24.94 ± 2.49 

29.71 ±2.97 

from Radon 

36.1 ±3.6 

34.0 ±3.4 

70.0 ±7.0 

Internal 

2.34 ±0.23 

13.83 ± 1.38 

16.17 ± 1.62 

2°®T1 Internal 

0.49 ±0.05 

2.93 ±0.29 

3.42 ±0.34 

i30Te 

0.12 ±0.02 

0.75 ±0.15 

0.87 ±0.17 

Total Expected 

43.8 ±3.7 

76.4 ±4.5 

120.2 ±8.1 

Data 

47 

76 

123 


TABLE VII: Numbers of expected background and observed 
two-electron events with Etot >2.0 MeV in Phases I and II, 
and for the combined data set, in the copper, natural tel¬ 
lurium, and ^^°Te foils. The combined data correspond to an 
exposure of 13.5 kg-yr. The contribution from 2z//3/3 decays 
of ^^°Te is negligible. 


VI. SEARCH FOR NEUTRINOLESS DOUBLE p 
DECAY 

The search for 0ul3f3 decays is performed by first se¬ 
lecting two-electron events using the criteria described in 
Section |IV[ where we require two electrons emitted from 
a common vertex in one of the ^^^Mo foils with a com¬ 
bined energy F^tot > 2 MeV. We then search for an excess 
in data above the background expectation in the F^tot 
distribution for energies close to the value of Q/ 3 / 3 . The 
contributions of the background from external sources, 
from radon, and from the internal ^^^Bi and ^^^Tl foil 
contaminations are fixed to the measured values given in 
Section |Vl 

We obtain the 2iyf3^ background contribution by fit¬ 
ting the F^tot distribution in the range F^tot > 2 MeV 
using the shape of the spectrum predicted by the Single 
State Dominance model for the 2^/3/3 decay of ^^^Mo [ 12 ]. 
The other background components are also taken into 
account in the fit. Figure shows that the fitted F^tot 
distributions for Phase I, Phase II, and for the combined 
data set agree with the data. The fitted number of 2^/3^ 
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FIG. 17: Distribution of the lengths of delayed a tracks for composite ^°°Mo foils for Phase II: (a,c) for electron and a, tracks 
on the same side of the foils, (b,d) for electron and a tracks on opposite sides of the foils, (a,b) for a tracks on the inner 
side of the foils, (c,d) for a tracks on the outer side of the foils. The data are compared to the simulated background with a 
normalisation determined by the ht of the different components of ^^^Bi background. “SW” corresponds to the ^^^Bi deposition 
on the surface of the wires, and “SF” to the deposition on the surface on the foil, where “IN” and “OUT” corresponds to the 
components from the wires and surfaces inside and outside relative to the position of the foil. “Internal” ^^^Bi contamination 
originates inside the ^°°Mo foils, and the “mylar” contamination from inside the mylar. 


events for E’tot > 2 MeV corresponds to a ^^^Mo half-life 
of 

Ti/2(2z// 3/3) = [6.93 ± 0.04 (stat)] x 10^^ yr, (9) 

after correcting for the signal efficiency, which is in agree¬ 
ment with the previously published result for Phase I [6] 
and with the world average [13]. 

The F^tot distribution in the region 2.8 < F^tot ^ 
3.2 MeV is shown in Figure and the different com¬ 
ponents of background in this energy window, and the 
number of observed two-electron events are given in Ta¬ 


ble VIII In Phase II, the observed background rate for 

2.8 < F^tot < 3.2 MeV is 0.44 ± 0.13 counts/yr/kg, with 
about 55% originating from 21 /decays of ^^^Mo, about 
20% from the radon gas contamination inside the track¬ 
ing chamber, and about 20% from internal ^^^Tl contami¬ 
nation in the ^^^Mo foils. We estimate the internal ^^^Bi 
contamination in the composite ^^^Mo foils to be 5%, 
while this background is negligible for metallic foils. The 
contributions from external backgrounds are also negli¬ 
gible. 

Since we observe no significant excess in data above 
the background expectation, a limit on the Oz//3/3 de¬ 
cay of ^^^Mo is derived. The uncertainties on the effi- 


Data Set Phase I Phase II Combined 


External Background 
^^^Bi from Radon 
2^^Bi Internal 
2°®T1 Internal 
21/(313 Decays 
Total Expected 
Data 


< 0.04 < 0.16 < 0.2 

2.8 ± 0.3 2.5 ± 0.2 5.2 ± 0.5 

0.20 ± 0.02 0.80 ± 0.08 1.0 ± 0.1 

0.65 ± 0.05 2.7 ± 0.2 3.3 ± 0.3 

1.28 ± 0.02 7.16 ± 0.05 8.45 ± 0.05 

4.9 ± 0.3 13.1 ± 0.3 18.0 ± 0.6 

3 12 15 


TABLE VIII: Numbers of expected background and observed 
two-electron events in Phases I and II in the ^°°Mo foil for an 
exposure of 34.3 kg-yr in the range Etot = [2.8 — 3.2] MeV. 
The Oz^/3/3 signal detection efficiency is 4.7% in this energy 
range. 


ciency to detect 0t/I3I3 events and on the estimated back¬ 
ground contributions are the two main components of 
the systematic uncertainty. As discussed in Section HYi 
the systematic uncertainty on the detection effi¬ 

ciency is 5%. The systematic uncertainties on the esti¬ 
mated background contributions are due to the activities 
of 2i/(3P decays, and the ^^^Bi and ^^^Tl backgrounds. 
An uncertainty of 0.7% on the 2i/(3P activity is obtained 
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FIG. 18: Distribution of the delayed a track length for metallic ^°°Mo foils (see Figure 17 caption for further details). 


from the fit to two-electron events in the energy range 
F^tot > 2 MeV. As discussed in Section |V| the systematic 
uncertainty on the normalisations of the background con¬ 
tributions from radon, and ^^^Tl radioactive con¬ 

taminants is 10%. This systematic uncertainty is taken 
into account in setting the limit on the decay of the 
lo^Mo isotope. The contributions of the external back¬ 
grounds and from thoron are negligible. 

The limit on the half-life is set using a modified 

frequentist analysis that employs a log-likelihood ratio 
test statistics m- The method uses the full informa¬ 
tion of the binned energy sum distribution in the F^tot = 
[2.0 — 3.2] MeV energy range for signal and background 
(see Figure as well as the statistical and systematic 
uncertainties and their correlations, and is described in 
more detail in mus]. All limits are given at the 90% 
C.L. The data are described well by the background-only 
hypothesis with a p value of p = 1 — CL}) = 0.647. Tak¬ 
ing into account the detection efficiency of 11.3% 

for the combined data set and the total exposure of 
34.3 kg-yr, we obtain a limit of Ti/2(Oz//3/3) > 1.1 xlO^"^ yr 
for the 0ul3f3 decays of ^^^Mo with decay kinematics sim¬ 
ilar to that for the light Major ana neutrino exchange. 

The result agrees with the median expected sensitivity 
of the experiment of Ti/2(0z^/5/3) = 1.0 x 10^^ yr within 
the ±1 standard deviation (SD) range of [0.7,1.4] x 
10^^ yr. This result is a factor of two more stringent 
than the previous best limit for this isotope [6] . The cor¬ 
responding upper limit on the effective Majorana neu¬ 


trino mass is (m^y) < 0.33-0.62 eV, where the range is 
determined by existing uncertainties on the calculations 
of the NMEs [161 ITStiST] and phase space factors [22123] . 
The upper value 0.62 eV is lower than the upper value 
previously reported in our rapid communication [7], be¬ 
cause of the use of the new NME calculation from [TB] . 
which is an update of the previous calculation m- 

We also derive constraints on other lepton-number vio¬ 
lating models: the supersymmetric models, the right-left 
symmetric models, and Majoron emission. 

In supersymmetric models, the process can be 

mediated by the exchange of a gluino or neutralino. Us¬ 
ing the obtained limit of Ti/2(0z//3;5) > 1.1 x 10^^ yr and 
the NME from [28] an upper bound is obtained on the 
trilinear i7-parity violating supersymmetric coupling of 
Am < (4.4 — 6.0) x 10“^/, where 


/ = 


Mq 

1 TeV 


Ms . 

1 TeV/ 


1/2 


( 10 ) 


and Mq and Mg represent the squark and gluino masses. 

Right-left symmetric models include right-handed cur¬ 
rents in the electroweak Lagrangian that predict different 
angular and energy distributions of the final state elec¬ 
trons from the decays. The NEMO-3 experiment, 

with the topological information for the two final-state 
electrons, can discriminate between the topologies from 
different mechanisms [24]. The corresponding half-life 
limits are given in Table [iX] and translate into an upper 
bound on the coupling between right-handed quark and 
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FIG. 19: Distribution of F^tot for two-electron events with F^tot > 2 MeV for the copper, ^^°Te, and natural tellurium foils 
(a,c,e), and for ^°°Mo foils (b,d,f), for Phase I (a,b) and Phase II (c,d), and combined (e,f). The combined data correspond to 
an exposure of 13.5 kg-yr for the copper, ^^°Te, and natural tellurium foils, and 34.3 kg-yr for the ^°°Mo foils. The data are 
compared to the sum of the expected background from 2i//3/3 decays of ^°°Mo, radon, external backgrounds, and from internal 
^^^Bi and ^°®T1 contaminations inside the foils. Only the 2z7/3/3 background contribution is fitted to the data, while the other 
background components are set to the measured values given in Section ^ Lower panels show residuals between data and 
expected background, normalized to the Poisson error, ignoring bins with 0 events. 
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lepton currents of (A) < (0.9 —1.3) x 10“^ and into an up¬ 
per bound on the coupling between right-handed quark 
and left-handed lepton currents of ( 77 ) < (0.5—0.8) x 10“^. 
The constraints are obtained using the NME calculations 
from [25H2Z]. 

The OuPP decay could also be accompanied by a Ma- 
joron (M), which is a light or massless boson that weakly 
couples to the neutrino [29]. In this case the energy 
sum of the two emitted electrons, Etot, will have a broad 
spectrum in the range [ 0 -( 5 / 3 / 3 ]. The shape will depend 
on the spectral index n, which determines the phase 
space dependence on the energy released in the decay, 
oc (Q/ 3/3 — The lower bound on the half- 

life of the OuPP decay with the spectral index n = 1 is 
given in Table [TX| The limit is set using the same method 
as the one used to extract the limit on the 0uf3l3 half- 
life with the energy sum of the two emitted electrons, 
Etot. This limit is almost a factor of two more strin¬ 
gent than the previous best limit for this isotope [30] . 
Taking into account the phase space factors given in m 
and the NME calculated in niiTHHii], an upper bound 
on the Majoron-neutrino coupling constant is obtained, 
(gee) < (1.6 - 3.0) X 10-5. 

The limits on lepton number violating parameters ob¬ 
tained here have comparable sensitivity to the best cur¬ 
rent results obtained with other isotopes, as shown in 
Table and in Eigure for the light Major ana neu¬ 
trino mass mechanism. 



Statistical 

Including Systematics 




Expected 

Oh'PP Mechanism 

Obs. 

Obs. 

-1 SD Median +1 SD 

Mass Mechanism 

1.1 

1.1 

0.7 1.0 1.4 

RH Gurrent (A) 

0.7 

0.6 

0.4 0.5 0.8 

RH Gurrent ( 77 ) 

1.0 

1.0 

0.6 0.9 1.3 

Majoron 

0.050 

0.044 

0.027 0.039 0.059 


TABLE IX: Observed and median expected lower limits on 
half-lives of lepton number violating processes (in units of 
10^^ yr) at the 90% C.L. using statistical and systematical 
uncertainties. The observed lower limits are also given using 
only statistical uncertainties. 


VII. CONCLUSIONS 

We have presented results based on an analysis of the 
full NEMO-3 data set with an exposure of 34.3 kg-yr of 
lo^Mo, which corresponds to 4.96 effective years of data 
collection and 6.914 kg of ^^^Mo. The calibration of the 
calorimeter, the long-term stability of data taking, and 
the determination of the backgrounds are discussed in 
detail. No evidence for 0uf3f3 decays of ^^^Mo has been 
found, as previously reported in our rapid communica¬ 
tion [7]. Taking into account statistical and systematic 
uncertainties, the limit on the 0ul3f3 decay half-life with 
decay kinematics similar to that for light Majorana neu¬ 
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FIG. 20: The 90% C.L. lower limits on Ti/ 2 ( 0 z^yd/ 3 ) for the 
light Majorana neutrino mass mechanism and upper limits 
on the effective Majorana neutrino mass {rrijy) using the same 
NME calculations [T^ [TsHST] and recent phase space calcula¬ 
tions [ 22 I [23] . The shaded regions correspond to the ranges 
from using different NME calculations. The hatched area 
corresponds to the expected range for {rriu), calculated from 
the neutrino oscillation parameters and assuming the inverted 
neutrino mass hierarchy. 


trino exchange is Ti/ 2 (Oz// 3 / 3 ) > 1.1 x 10^^ yr (90% C.L.). 
The corresponding limit on the effective Majorana neu¬ 
trino mass is in the range (m^y) < 0.33-0.62 eV, depend¬ 
ing on the NME calculation used in the derivation. 

Studies of the backgrounds using various decay chan¬ 
nels, radioactive sources, and HPGe measurements be¬ 
fore the installation of the detector are used to con¬ 
struct and validate a detailed model of the background. 
In Phase II, the expected background rate in the 
signal region E^tot = [2.8 — 3.2] MeV is 0.44 ± 
0.13 counts/yr/kg. About half of this background is ex¬ 
pected to be 2uPf3 decays of ^^^Mo, and the remaining 
background is caused in roughly equal parts by the radon 
gas contamination inside the tracking chamber, which is 
about 5 mBq/m^, and by ^^^Tl contamination inside the 
ioomo foils, which is between 90-130 /iBq/kg depending 
on the type of foil. No background events are observed in 
the region of Etot = [3.2 — 10] MeV for NEMO-3 sources 
containing isotopes with < 3.2 MeV (^^^Mo, ^^Se, 
^^^Te, ^^^Cd), or in the copper foil, which is not a double 
(3 emitter, during the entire running period correspond¬ 
ing to an exposure of 47 kg-yr. 

This low level of background demonstrates that an ex¬ 
tremely low level of non double P background can be 
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Half-Life 
( 10 ^® yr) 

(eV) 

Yec 

(eV) 

(A) 

( 10 -®) 

{n) A'lii/f 

( 10 -®) ( 10 - 2 ) 

{ 9 ee) 

( 10 -®) 

[This Work] 

0.11 

0.33-0.62 

0.33-0.62 

0.9-1.3^ 

0.5-0.8^ 4.4-6.0 

1.6-3.0^ 

[3l|33] 

0.28 

0.3-0.71 

0.31-0.75 

1.6-2.4^ 

0.9-5.3^ 

17-33^ 

[MES] 

1.9 

0.14-0.34 0.14-0.34 



0 . 8 - 1.6 

’'®Ge [36] 

2.1 

0.2-0.4 

0.26-0.62 




’'®Ge |37l|38] 

1.9 

0.35 

0.27-0.65 

1.1 

0.64 

8.1 


TABLE X: Limits at the 90% C.L. on half-lives and lepton number violating parameters. Published experimental constraints 
on {rrijy) and recalculated values with NMEs from Refs. [161 [TSHIII |39] 9 ^ 1 so given. 


obtained with half-lives in Table 


IX 


using the half-life limit of 2.1 x 10 ^^ yr, using the half-life limit of 2.2 x 10 ^^ yr. 


achieved by the future SuperNEMO experiment, which 
will employ the NEMO-3 technique. The SuperNEMO 
Collaboration proposes to search for 0ul3f3 decays using 
100 kg of double P isotopes [24]. The 21 /background 
will be further reduced by improving the energy reso¬ 
lution and by measuring an isotope with a long 2i/PP 
half-life, currently assumed to be ^^Se. Other favorable 
isotopes, such as ^^^Nd and ^^Ca, are also studied. A 
first SuperNEMO demonstrator module, currently under 
construction, will contain 7 kg of ^^Se. The objective 
is to demonstrate that the background can be reduced 
by 1-2 orders of magnitude compared to the NEMO-3 


detector. 


Acknowledgments 

The authors would like to thank the Modane Under¬ 
ground Laboratory staff for their technical assistance in 
running the experiment. We acknowledge support by the 
grants agencies of the Czech Republic, CNRS/IN2P3 in 
France, RFBR in Russia, STFC in U.K. and NSF in U.S. 


[1] P. Minkowski, Phys. Lett. B 67, 421 (1977); T. Yanagida, 
Prog. Theor. Phys. 64, 1103 (1980); M. Gell-Mann, P. 
Ramond and R. Slansky, Proceedings of the Supergrav¬ 
ity Stony Brook Workshop, New York 1979, eds. P. Van 
Nieuwenhuizen and D. Freedman; R. N. Mohapatra and 
G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[2] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 
(1986). 

[3] J. Schechter and J.W.F. Valle, Phys. Rev. D 25, 2951 
(1982). 

[4] R. Arnold et al, Nucl. Inst. Meth. A 536, 79 (2005). 

[5] S. Rahaman et al, Phys. Lett. B 662, 111 (2008). 

[ 6 ] R. Arnold et a/., Phys. Rev. Lett. 95, 182302 (2005). 

[7] R. Arnold et al, Phys. Rev. D 89, 111101 (2014). 

[ 8 ] GEANT Detector Description and Simulation Tool, GERN 
Program Library Long Writeup W5013 (1995). 

[9] O.A. Ponkratenko, V.I. Tretyak, and Yu.G. Zdesenko, 
Phys. Atom. Nucl. 63, 1282 (2000). 

[10] J. Argyriades et al., Nucl. Inst. Methods A 606, 449 
(2009). 

[11] R. Arnold et al., Phys. Rev. Lett. 107, 062504 (2011). 

[12] J. Abad et al., Ann. Fis. A 80, 9 (1984); P. Domin et al., 
Nucl. Phys. A 753, 337 (2005). 

[13] Particle Data Group, K.A. Olive et al., Ghin. Phys. G 
38, 090001 (2014). 

[14] T. Junk, Nucl. Inst. Methods A 434, 435 (1999). 

[15] J. Argyriades et al., Phys. Rev. G 80, 032501 (2009). 

[16] J. Hyvarinen and J. Suhonen, Phys. Rev. G 91, 024613 
(2015). 

[17] J. Suhonen and O. Givitarese, J. Phys. G 39, 124005 

( 2012 ). 

[18] F. Simkovic, V. Rodin, A. Faessler, and P. Vogel, Phys. 


Rev. G 87, 045501 (2013). 

[19] J. Barea, J. Kotila, and F. lachello, Phys. Rev. G 91, 
034304 (2015). 

[20] P.K. Rath et al., Phys. Rev. G 88 , 064322 (2013). 

[21] T.R. Rodriguez and G. Martinez-Pinedo, Phys. Rev. 
Lett. 105, 252503 (2010). 

[22] J. Kotila and F. lachello, Phys. Rev. G 85, 034316 (2012). 

[23] S. Stoica and M. Mirea, Phys. Rev. G 88 , 037303 (2013); 
Updated in arXiv: 1411.5506 

[24] R. Arnold et al., Eur. Phys. J. G 70, 927 (2010). 

[25] J. Suhonen, Nucl. Phys. A 700, 649 ( 2002 ). 

[26] T. Tomoda, Rept. Prog. Phys. 54, 53 (1991). 

[27] K. Muto, E. Bender, and H.V. Klapdor, Z. Phys. A 334, 
187 (1989). 

[28] A. Faessler, G.L. Fogli, E. Lisi, A.M. Rotunno, and F. 
Simkovic, Phys. Rev. D 83, 113015 (2011). 

[29] P. Bamert, G.P. Burgess, and R.N. Mohapatra, Nucl. 
Phys. B 449, 25 (1995). 

[30] R. Arnold et al., Nucl. Phys. A 765, 483 (2006). 

[31] M. Doi, T. Kotani, and E. Takasugi, Phys. Rev. D 37, 
2575 (1988). 

[32] E. Andreotti et al., Astropart. Phys. 34, 822 (2011). 

[33] G. Arnaboldi et al, Phys. Lett. B 557, 167 (2003). 

[34] A. Gando et al., Phys. Rev. Lett. 110, 062502 (2013). 

[35] A. Gando et al., Phys. Rev. G 86 , 021601 ( 2012 ). 

[36] M. Agostini et al., Phys. Rev. Lett. Ill, 122503 (2013). 

[37] H.V. Klapdor-Kleingrothaus et al., Eur. Phys. J. A 12 , 
147 (2001). 

[38] M. Gunther et al., Phys. Rev. D 55, 54 (1997). 

[39] J. Menendez, A. Poves, E. Gaurier, and F. Nowacki, Nucl. 
Phys. A 818, 139 (2009). 









